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ABSTRACT 

We analyze interferometric measurements of the Luminous Blue Variable Eta Carinae with the 
goal of constraining the rotational velocity of the primary star and probing the influence of the 
companion. Using 2-D radiative transfer models of latitude-dependent stellar winds, we find that 
prolate wind models with a ratio of the rotational velocity (v rot ) to the critical velocity (v CT [ t ) of 
W = 0.77 - 0.92, inclination angle of i = 60° - 90°, and position angle PA = 108° - 142° reproduce 
simultaneously K-band continuum visibilities from VLTI/ VINCI and closure phase measurements 
from VLTI/ AMBER. Interestingly, oblate models with W = 0.73 - 0.90 and i = 80° - 90° produce 
similar fits to the interferometric data, but require PA = 210° — 230°. Therefore, both prolate and 
oblate models suggest that the rotation axis of the primary star is not aligned with the Homunculus 
polar axis. We also compute radiative transfer models of the primary star allowing for the presence 
of a cavity and dense wind- wind interaction region created by the companion star. We find that 
the wind-wind interaction has a significant effect on the if -band image mainly via free-free emission 
from the compressed walls and, for reasonable model parameters, can reproduce the VLTI/ VINCI 
visibilities taken at </> v t>03 = 0.92 — 0.93. We conclude that the density structure of the primary 
wind can be sufficiently disturbed by the companion, thus mimicking the effects of fast rotation in 
the interferometric observables. Therefore, fast rotation may not be the only explanation for the 
interferometric observations. Intense temporal monitoring and 3-D modeling are needed to resolve 
these issues. 

Subject headings: stars: atmospheres — stars: mass-loss — stars: variables: general — supergiants 
— stars: individual (Eta Carinae) — stars: rotation 



1. INTRODUCTION 

Eta Carinae is one of the most luminous objects in the 
Galaxy, allowing for the study of massive stellar evolu- 
tion under extreme conditions. Eta Car also presents us 
with a unique opportunity to witness the early evolution 
of violent giant outbursts that happen during the Lu- 
minous Blue Variable (LBV) phase, like that which cre- 
ated Eta Car's massive H omunculus nebula (jThackeravl 
119501 : ISmith et al1l2003bl). Eta Car is general ly accepted 
to be a binary system (jDamineli et al.|["l997[ ) comprised 
of two massive stars, t]a (primary) and t]b (secondary). 
The total luminosit y is dominated by 77 a and amou nts to 
Ltot > 5 x 10 6 L Q ([Davidson fe Humphreys! 1 1 99 7T). and 
the total mass is of at least 110 M (jHillier et al.l [2001. 
hereafter H01). Most authors a gree on a high eccentricity 
(e ~ 0.9, iCorcoran et al.ll2001[) and an orbital period of 
2022.7 ± 1.3 davs (jDamineli et al.ll2008D . while the other 
orbital parameters are uncertain. 

t]a is an LBV star with a high mass-loss rate of 
M ~ 1O _3 M0 yr _1 and wind terminal velocity of 
Voo ~ 500 - 600 km s" 1 (|H01h . The stellar parame- 
ters and evolutionary state of t]b are uncertain since it 
has never been directly detected, and only indirect in- 
formation is available. A broad range of t]b stellar pa- 
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rameters can explai n the ionization of the inner eject a 
(jMehner et al.|[20Tol ). and X-ray observations require tjb 
to ha ve ~ 3000 km s" 1 a nd M ~ 10~ 5 M© yr" 1 
(e.g.. IPittard fc Corco ran 2002). 3-D numerical simula- 
tions show that the wind of t]b influences the geometry 
of the very dense wind of t]a via th e creation of a cavity 
and dense wind- wind collision zone (IPittard fc Corcoran! 
2002; lOkazaki et "all 120081 : iParkin et al.ll2009D . 

Based on the variations of Ha absorption line pro- 
files in scattered light from the Homunculus, which pro- 
vide us with different viewing directions to the star, 
and assuming that the rotati on axis of 77 a and t he Ho- 
munculus axis are aligned, ISmith et al.l (|2003ai here- 
after S03) suggested that t]a has a latitude-dependent 
wind, with faster, denser outflow in polar directions. In- 
terferometric measurem ents obtained wit h VLTI /VINCI 
(jvan Boekel et al.l 120031, hereafter VB03 ; iKervellal l2QQ7h 
and VLTI/ AMBER ([Weigelt et alJl2QQ7l hereafter W07) 
in the near-infrared if -band continuum are consistent 
with a n e llipsoi dal shape projected on the sky. Both 
VB03| and IW07I interpreted this as evidence for a dense 
prolate wind generated by f ast ro t ation, as theoretically 
predicted by Ow ocki et al.l (|1996l Il998h . However, the 
influence of t]b on the Ha absorption profiles and on the 
if -band emission, which are the key diagnostics support- 
ing the fast rotation of t?a, is poorly constrained. 

We analyze published K-band continuum interferomet- 
ric data of Eta Car obtained at orbital phased v bO3 = 
0.92-0.93 (fVB03h and wO7 = 0.27-0.30 (|W07h . Our 
goals are to constrain the rotational velocity and spatial 

5 Assuming the ephemeris from Damineli et al. (2008). 
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orientation of the rotation axis of tja based on the effects 
of rotation on the wind density structure, which deter- 
mines the geometry of the if -band emitting region. We 
also investigate the influence of tjb on the inner wind of 
t]a through the presence of a low-density cavity and den- 
sity enhanced wind-wind collision zone, with the goal of 
determining how these may affect the interpretation of 
the VINCI dataset obtained relatively close to periastron. 
We apply, for the first time for Eta Car, 2-D radiative 
transfer models of latitu de-de pendent winds generated 
by rapid rotation (Section l2TTj) and of a modified wind of 
r]A which includes a w ind cavity and coll idin g wi nd inter- 
action region (Section l3.ip . Sections 12 . 21 and [3721 compare 
the computed interferometric observables for each sce- 
nario with the available observations. Section [4] discusses 
whether t]a is a fast rotator, a possible misalignment be- 
tween its rotation axis and the Homunculus polar axis, 
and the influence of t]b on the wind of t?a- 

2. SINGLE-STAR, LATITUDE-DEPENDENT WIND 
GENERATED BY RAPID ROTATION 

2.1. Radiative transfer modeling 

We compute 2-D latitude-dependent wind models for 
r]A using an updated vers ion of the 2-D radiative transfer 
code of iBusche fc Hillierl (|2005l hereafter BH 05) and via 
a met hodology similar to that described in iGroh et"aT] 
(2008). In the following, we briefly describe only the 
essent i al aspects of the code. We r efer t he reader to 
BH05L IGroh et al.l (|200d I20Q8L I2009ah , and lDriebe et al.l 
(|2009[ ) for further details on the code's applications. 

The 2-D models use as input several quantities 
(e.g., energy- level populations, ionization structure, 

and radiatio n fiel d) from the spherically symmetric 

model of t]a (HOlj; lHQ6h computed using the non-LTE, 
fully line blanke ted radiative transfer code CMFGEN 
(Hillier & Miller 1998). We assume the same parame- 
ters derived by stellar temperature T* = 35,310 K 
(at Rosseland optical depth tr oss = 150), effective tem- 
perature T e ff= 9,210 K (at tross = 2/3), luminosity 
L* = 5xl0 6 L©, M = 10~ 3 M© yr-\ Voo = 500km s _1 , 
a clumping volume-filling factor / = 0.1, and distance 
d = 2.3 kpc. The a tomic model and abundances are 
described by EoT| and IH06I . 

The code allows for the specification of any arbitrary 
latitude-dependent variation of the wind density p and 
Voo. For the l atitudinal vari a tion o f p, we adopt the pre- 
dictions from lOwocki et al.l (jl998) for gravity-darkened 
(GD) line-driven prolate winds, 
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and the non-gravity-darkened (non-GD) scaling for 
oblate winds, 
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where is the colatitude angle (0°= pole, 90°= equa- 
tor), po is the density at the pole, and W is the ratio 
of the rotational velocity of t]a (v ro t) to its critical ve- 
locity (^crit)- Note that t he oblate-wind density scaling 
assumes a standard CAK ([Castor et al.lll975[ ) a param- 
eter value of 2/3. The inclination angle i is the tilt an- 
gle between the rotation axis of the star and the line 



of sight (i = 0°= pole-on view, i = 90°= equator-on 
view). Emissivities and opacities are scaled according to 
the modified density, but the code assumes a spherically 
symmetric ionization structure. Gravity darkening and 
distortion of the shape of the hydrostatic core of the star 
are not included at this time; however, these are not ex- 
pected to significantly affect the if -band emitting region 
since the source function in the infrared continuum de- 
pends only weakly on the temperature. For consistency, 
we also assume the lOwocki et al.l (|1998[ ) predictions for 

Voo(0). 

We calculated a grid of prolate and oblate models with 
i ranging from 0° to 90° in steps of 1°, and W ranging 
from to 0.99 in steps of 0.01, totaling 18,200 models. 
For a desired position angle (PA) orientation on the sk}0, 
the if -band image, visibilities, and closure phases (CP) 
are computed. 

2.2. Results 



Since the spherical CMFGEN model of H01 d oes no t 
fit the geometry of the if -band continuum (VB03; W07), 
we first attempted prolate wind models (Figure [1]) to re- 
produce the suggested elongation of the if -band contin- 
uum observed with VINCI in 2003 January-February. 
We find the model parameters to be degenerate, as mod- 
els with PA = 108°-142°, W = 0.80-0.99, and i = 30°- 
90° fit the VINCI data equally well. Figure [2^ presents 
a color-coded plot of the reduced x 2 values of the fit to 
the VINCI data, with the whitish regions corresponding 
to the best-fit models. In general, models with relatively 
large W require lower z, while those with relatively lower 
W need larger i to fit the data. 

Additional constraints to W and i can be obtained by 
analyzing the closure phase, which mea sures the asym- 
metry of the brightness distribution. IW07I measured 
CP= 0° ± 3° in the K-band continuum in 2005 Feb, and 
a color-coded plot of the model CPs, computed using 
the same telescope array configuration used by IW071 is 
shown in Figure [2)3 . We find that prolate models with 
high W (0.82 to 0.99) and intermediate inclination an- 
gles (i ~ 30° — 60°) yield a K-band image which is no- 
ticeably not point-symmetric (Figure Q}i) , and thus are 
ruled out since they have much larger CPs (7° to 19°) 
than observed. To fit simultaneously the VLTI/ VINCI 
visibilities and the VLTI/ AMBER CPs, W = 0.77-0.92 
and i = 60° — 90° are needed (Figure [2k). Strikingly, 
this inclination is signific antly higher than that of the 
Homunculus (iuom =41°; [Smith! 12006), which has been 
assumed so far to align with the current rotation axis of 
t]a- We find that, based on the available interferometric 
data, if t]a has a prolate wind, its rotation axis is not 
necessarily aligned with the Homunculus polar axis. 

The image asymmetries from models with i = 30° 
to 60° (red spot in Figure [TJi) arise from the latitude- 
dependent K-band photosphere associated with the wind 
density variations. Rays from the NW hemisphere cross 
the low-density equator, and so accumulate a lower opti- 
cal depth than those from the SE hemisphere at the same 
impact parameter. The deeper penetration of NW vs. SE 
rays thus leads, for the nearly spherical source function 
which is larger for smaller radii, to a stronger emerging 

6 Measured in degrees east of north. 
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Figure 1. Left panels: Density structure of the latitude-dependent wind models of t]a in the x-z plane (i. e., equator-on). Middle: K-band 
image projected on the sky. Right: VINCI visibilities for the 24 m baseline as a function of telescope PA (connected black asterisks) 
compared to the respective model prediction (red dotted line). Note that the projected baseline length of the VINCI measurements changes 
as a function of PA. From top to bottom, the panels correspond to examples of best-fit prolate wind models with W = 0.85, i = 75°, 
PA = 130° (panels a — c) and W = 0.92, i = 80°, PA = 108° (panels d — f), and prolate wind aligned with the Homunculus but with large 
closure phase (panels g - i; W = 0.99, i = 41°, PA = 130°). 



radiation field, and thus a higher surface brightness in 
the NW hemisphere. 

One may ask if a prolate wind model is unique in being 
able to fit the interferometric data. Therefore, we also 
analyze the visibility as a function of telescope PA and 
CP for oblate wind models. We find that the oblate 
model fits the visibilities as well as the prolate model, 
provided that PA = 210° - 230° (Figure [3]), i.e., the 
stellar rotation axis projected on the plane of the sky 
is rotated by 80°-100° from the Homunculus symmetry 
axis. Oblate models with W = 0.73-0.90 and i = 80 - 
90° are able to fit both the observed visibilities and CP 
(Figure [2jl-f). Therefore, for oblate wind models, the 
rotation axis of t]a is not aligned with the Homunculus 
polar axis. 

3. EFFECTS OF THE WIND- WIND COLLISION ZONE 
CREATED BY THE COMPANION 

Since Eta Car is believed to be a colliding-wind binary 
system, we also investigate the effect of the low-density 
cavity and dense wind- wind interaction region created by 
the wind of t]b in the wind of t]a- 



3.1. Radiative transfer modeling 

Our models are based on the spherically symmetric 
models of t]a (HO I; H 06|), with the same stellar par am e- 
ters as in Section |2TT| but using the 2-D code of BHOl to 
create a low-density cavity and dense interaction-region 
walls. Further details are given in Groh et al. 2010 (in 
preparation), and here we outline only the main charac- 
teristics of our implementation. 

We approximate the cavity as a conical surface with 
half-opening angle a and interior density 0.0016 times 
lower than that of the spherical wind model of t^a- We 
include cone walls of angular thickness 5a and, assum- 
ing mass conservation, a densit y contrast in th e wall 
of f a = [1 - cos(a)]/[sin(oOH (|Gull et all [2009( 1 times 
higher than the wind density of the spherical model of 
r]A at a given radius. The conical shape is justified since 
the interferometric observations were taken at orbital 
phases sufficiently before periastron when such a cav- 
i ty has an approxima tely 2-D axisymmetric conical form 
(Oka zaki et al.l 12008). Based on the expected location 
of the cone apex during these phases ( Can to et al.l ll996; 
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Figure 2. Left column: From top to bottom, color-coded plots of 
the reduced x 2 values of the fit of the VINCI visibilities (a), closure 
phase (b), and reduced % 2 values rejecting models with CP > 5° 
(c), as a function of W and i, for prolate models with PA = 130°. 
Right column: Idem, for oblate models with PA = 220° (d-f). 

lOkazaki et al.l [2QQ8h . we place the cavity at a distance 
^apex from the primary star. We assume that the mate- 
rial inside the cavity and along the walls has the same 
ionization structure as the wind of t?a, implying that the 
shocked wind of the primary star cools radiatively. 

3.2. Results 

The amount of influence of the wind- wind interaction 
on the observables depends on two factors: how close the 
cavity gets t o the K-band emitting re gion of t]a ("bore 
hole effect", Ma dura fc OwocO 12010). and how much 
free-free radiation is emitted by the dense walls of the 
shock cone ("wall effect"). We find that the interfero- 
metric observables in the K-band are insensitive to the 
bore hole effect if <i a pex ^ 8 AU. This is because most of 
the observed K-band emission comes from a region with 
a cha racteristic 50% encircled-energy radius of 4.8 AU 
(|W07l ). and a halo with significant emission exists up to 
~ 8 AU. For the st andard assumed orb ital parameters 
of Eta Car (see, e.g.. lOkazaki et aLll2008l ). <i a pex ^ 8 AU 
during most of the orbit (0.055 < (j) < 0.945), which 
includes the epochs when the VINCI (</> v b03 = 0.92- 
0.93) and AMBER (0 wO7 = 0.27-0.30) observations were 
taken. The VINCI observations were obtained when 



<4 



10 AU, thus a bore hole effect alone is not able 



to explain the VINCI observations. 

According to our models, the main effect from the 
wind-wind collision zone at the orbital phases analyzed 
is extended free- free emission from the compressed walls. 
Such an effect may influence the geometry of the K-band 
emitting region even at orbital phases far from perias- 
tron, depending on the observer's location, geometry of 
the cavity, z, a, a nd f n . A model w ith a = 54°, f a = 9.7 
(i. e., 5a = 3°; IGull et aLll2QQ9h . and i = 41° shows 
significant wall emission and is able to produce a signifi- 
cant elongation of the K-band continuum (Figure HJ. To 
explain the VINCI observations, the symmetry axis of 
the cavity has to be oriented along PA~ 35° — 45° (i.e., 
SW-NE axis; Fig. [4]), which is roughly consistent with 
that expected from a longitude of periastron of uo = 243° 
(e.g, Okazaki et al. 2001 iParkin et al.l 120091 IGull et al.l 
2009; Gro h et al.ll20Toh and a counterclockwise motion 
of t]b on the sky (Figure (Dd). The M of 77 a needs to be 
slightly reduced to 8 x 10 -4 M yr _1 to compensate for 
the extra extension in the K-band emitting region caused 
by the wall effect. Comparing Figures [T] and HJ the cavity 
and walls have an effect on the available interferometric 
observables that is as large as that due to a latitude- 
dependent wind caused by rapid rotation, although the 
morphology of the K-band images is noticeably different. 

The cavity model assumed here is very idealized, since 
it do es not take into acco unt the distortion of the shock 
cone (Okazaki et al. 2008), or the ins tabilities thought t o 
arise in the wind- wind collision zone (Parkin et al. 2009). 
Nevertheless, our conclusion is that, using a 2-D model, 
noticeable elongation of the K-band region can occur if 
the apex of the cone penetrates a significant distance into 
the K-band emitting region, that cone is open toward the 
observer, and/or the compressed walls are dense enough 
and emit free-free radiation. Three-dimensional radiative 
transfer models, when available, would be desirable to 
test the conclusions found here. 

4. DISCUSSION 

r]A is routinely referred to as the prototype of a massive 
star with a fast, dense polar wind created by rapid stellar 
rotation, even though the system is believed to contain 
a massive companion, t]b- How much does t]b affect the 
K-band emitting region of t/a? We show that, assuming 
the standard orbital and wind parameters of Eta Car, 
even if t]a has a spherical wind, its inner density struc- 
ture can be sufficiently disturbed by 7?b, mimicking the 
effects of a prolate/oblate latitude-dependent wind in the 
available interferometric observables in the K-band con- 
tinuum. Therefore, fast rotation may not be the only 
explanation for the interferometric observations. 

If we ignore the presence of the companion, our 2- 
D modeling shows that both single star prolate and 
oblate wind models are able to explain the elongation 
of the K-band emitting region of t]a and the CP mea- 
surements. While a prolate wind is thought to arise 
in a gravity- darkened, fast-rotat i ng; st ar with a radia- 
tive envelope (jOwocki et al.| [l996. 1998), an oblate wind 
can be produced by a fast -rotating star when gravity- 
darkening: is not imp ortant (IB i or kman fc Ca ssinelli 1993; 
lOwocki et all 119941 . 119961 . 119981 : llgnace et al.l Il996f ). for 
instance if the envelope is predominantly convective. 
Interestingly, convection can be present in the enve- 
lope of stars close to the Eddington limit (Langer 1997; 
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Figure 3. Similar to Figure^ but for a best-fit oblate wind model with W = 0.84, i = 85° , PA = 220°. 
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Figure 4. Panels a,c,d: Similar to Figure [T] but for the t]a model with a spherical wind including a cavity and compressed walls created 
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by the wind of tjb- The model is appropriate for the VINCI observations (</> = 0.93) and assumes i = 41°, PA 
a = 54°, and 5a = 3°. Panel b: Sketch of the binary orbit (yellow) and the shock cone orientation at <f> = 0.93 relative to the Homunculus 
(not to scale), assuming a counterclockwise motion of t]b on the sky. The orbital plane is assumed to be in the skirt plane with the orbital 
axis (blue) aligned to the Homunculus axis of symmetry. North is up and East is to the left. 



iMaeder et al.l l2QQ8h . such as t?a, although it is unclear 
whether convection is efficient enough to transport a sig- 
nificant fraction of the flux. 

We find that moderately fast rotation and high inclina- 
tion angles are required to fit simultaneously the VINCI 
and AMBER data, with W = 0.77-0.92 and i = 60° -90° 
for the best prolate models, and W = 0.73-0.90 and 
i = 80° — 90° for the best oblate models. The possi- 
ble tilting between the current rotation axis of t]a and 
the Homunculus axis is striking. If the Giant Eruption, 
which created not only the Homunculus ne bula, but als o 
a blast wave with much higher velocities (Smith 2008), 
happened in that would indicate that either the an- 
gular momentum loss was not axisymmetric, that strong 
dynamical interactions between t]a and t]b (or a third 
member) occurred in the mean time, or both. 

Alternatively, a speculative scenario that permits i = 
iHom =41° is one that allows for time variability in W, 
which means a different latitude dependency of the pro- 
late wind at the different epochs of the observations. This 
would require W = 0.93 - 0.99 during the VINCI obser- 
vations in 2003 January-February (0 v bO3 = 0.92-0.93), 
and a significantly decreased W = 0.75 — 0.82 during 
the AMBER observations in 2005 Feb (0 w q 7 = 0.27-0.30 
). A possible explanation for this variability could be a 
change in the stellar radius due to S-Dor type variability, 
which would cause v rot and W to vary, as has been seen 
for other LB Vs such as AG Car (jGroh et al.ll2006ll2009bD 
and HR Car (iGroh et al.ll2009afl 0. Regardless of such a 

7 A variability in v TO t, with a gradual increase before the spec - 
troscopic event, has been suggested before (S03; Davidson 1999), 



speculative scenario, based on the available inter feromet- 
ric data alone, we are unable to determine whether 77 a 
possesses an oblate or prolate wind, and whether or not 
the rotation axis is aligned with the Homunculus sym- 
metry axis (both i and PA). The fact that the wind- 
wind collision zone can mimic the effects of a latitude- 
dependent wind further complicates matters. Thus, more 
interferometric observations are needed to resolve these 
issues. 

Much of the support for the presence of a fast, dense 
polar outflow comes from spectroscopic observations of 
the Ha abso rption profile at diff erent locations in the 
Homunculus ( S03; iHillier fc Allerill992h . Assuming that 
the r otati on axis of t]a is aligned with the Homunculus 
axis, S03 found that the absorption profile of Ha ex- 
tends to velocities as high as —1200 km s _1 in the polar 
spectrum, while at mid-latitudes a value of — 600 km s _1 
is found. However, as pointed out by IS03L the latitu- 
dinal behavior of the Ha absorption can be affected by 
the presence of tjb if it is able to photoionize the region 
where H absorption would occur. Determining how much 
the wind terminal velocity changes as a function of lati- 
tude is only possible when the influence of t]b is properly 
determined. The exact influence of t]b is unknown and 
depends mainly on the number of its H ionizing pho- 
tons and the orbital parameters. It would be intrigu- 
ing to investigate how t]b may affect these Ha diagnos- 
tics. Detailed multi-dimensional radiative transfer mod- 
els of spectral lines and other diagnostics are urgently 
warranted, and will be the subject of future work. 

but in a different context (shell ejections). 
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